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RARE CHARM DECAYS 



EUGENE GOLOWICH 
Department of Physics, University of Massachusetts 
Amherst MA 01003, USA 



This paper is a written version of a talk on rare (FCNC) D meson decays as presented at 
the Electroweak Moriond 2002. The presentation proceeds in two parts. We first consider 
Standard Model predictions, taking into account both short-distance and long distance ef- 
fects. Then several New Physics options {e.g. supersymmatry, strong dynamics, extra large 
dimenions, etc) are considered. 



1 Introductory Remarks 

I dedicate this paper (as I did my talk) to Jean Tran Thanh Van and his wife for their valuable 
contributions over many years to our disciDline. 

My talk summarized a recent paper El written in collaboration with Gustavo Burdman, 
JoAnne Hewett and Sandip Pakvasa on the topic of rare flavor-changing neutral current (FCNC) 
decays of D mesons. Transitions of this class include: (i) D {V = p,uj, (p, . . .), ^ (ii) 

D ^ X£+r {X = 7r,K,r],p,L0,(P,...), (iii) D ^ Xu^u, {X = tt, K,i], . . .), (iv) D ^ 77, (v) 

D ^i+e-. 

Perhaps the most noteworthy aspect of the above list is that, at the time of this writing, 
exactly zero FCNC events have been detected. □ We are optimistic that this situation will be 
rectified in the fairly near future. The database for charm-related physics continues to expand, 
and experiments at B-factories, hadron collidiers and tau/charm facilities will soon achieve truly 
interesting sensitivities. 

As such, one can reasonably imagine a future conference talk which announces the detection 
of a D meson FCNC transition. How to interpret such a signal? There is a commonly adopted 
procedure which should be followed. First one determines whether the Standard Model (SM) 
can explain the observed events. If not, one proceeds to consider a menu of New Physics models. 



^The single-photon decays were considered previously by the authors in a separate paper. I 
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Figure 1: Dilepton mass distribution for D"*" ^ ir'^e^e . The daslied (solid) line is the short-distance (total) SM 
contribution. The dot-dashed line is the R-parity violating contribution. 



2 Standard Model Analysis 



In principle, the task of producing SM predictions for FCNC D meson decays is straightforward. 
There are two components to the analysis, short-distance (SD) and long-distance (LD), which 
must be separately calculated. We consider each in turn. 



2.1 The short- distance component 

Short-distance amplitudes are concerned with the QCD degrees of freedom (quarks, gluons) and 
any relevant additional fields (leptons, photons). Thus, the short distance part of the D — 
Xul'^(.~ amplitude involves the quark process c u£~^i~ . It is usually most natural to employ 
an effective description in which the weak hamiltonian is expressed in terms of local multiquark 
operators and Wilson coefficients. @ For example, the effective hamiltonian for c u£^£~ with 
renormalization scale fi in the range irib > fJ. > mc is^ 
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In the above, O^^l four-quark current-current operators, 0'^_q are the QCD penguin op- 
erators, O7 (Og) is the electromagnetic (chromomagnetic) dipole operator and Og^io explicitly 
couple quark and lepton currents. For example, we have 



2 

O; = ^m,(uLC7^,c^)F'^^ 0', = r^{uLl^CLm^e) . (2) 

The famous Inami-Lim functions! contribute to the Wilson coefficients C7-10 at scale fi = Mw- 
Figure ^ displays the predicted dilepton mass spectrum for — > ir^i^i^ . Several distinct 
kinds of contributions are included. The short-distance SM component corresponds to the dashed 
line, which is seen to lie beneath the other two curves. For reference, we cite the inclusive 'short 
distance' branching ratio, 

K'LxU^e- - 2 X 10-« . (3) 



''Quantities with primes have had the explicit 6-quark contributions integrated out 



(a) (b) 
Figure 2: Long distance contributions to — > i~ . 
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(a) (b) 
Figure 3: Long distance contributions to _D+ n^l^l^ . 

2.2 The long-distance component 

The long-distance component to a transition amplitude is often cast in terms of hadronic entities 
rather than the underlying quark and gluonic degrees of freedom. For charm decays, the long- 
distance amplitudes are typically important but difficult to determine with any rigor. There 
are generally several long-distance mechanisms for a given transition, e.g. as indicated for 
l+l- in Fig. | and for D+ tt+uu in Fig. |[ 
Let us return to the case of Tr^£'^£~ depicted in Fig. |l|. The solid curve represents 

the total SM signal, summed over both SD and LD contribution. In this case the LD component 
dominates, and from studying the dilepton mass distribution we can see what is happening. The 
peaks in the solid curve must correspond to intermediate resonances ((/>, etc). The corresponding 
Feynman graph would be analogous to that in Fig. |3|(a) in which the final state neutrino pair is 
replaced by a charged lepton pair. One finds numerically that 

R^(SM) _ m _ 6 /.N 



2.3 The Standard Model Predictions 

Basing our analysis in part on existing literature, i we have calculated both SD and LD ampli- 
tudes for a number of FCNC D transitions. El Results are collected in Table |l]. As stated earlier, 
the current database for processes appearing in Table |l] consists entirely of upper bounds (or in 
the case of — > 77 no data entry at all). In all cases existing experimental bounds lie below 
the SM predictions, so there is no conflict between the two. For some {e.g. D — > ■k£~^£~) the gap 
between SM theory and experiment is not so large and there is hope for detection in the near 
future. In others {e.g. £~^£~) the gap is enormous, leaving ample opportunity for signals 

from New Physics to appear. This point is sometimes not fully appreciated and thus warrants 
some emphasis. It is why, for example, attempts to detect AMd via D^-D^ mixing experiments 
are so important. 



Table 1: Standard Model predictions and current experimental limits for the branching fractions due to 
short and long distance contributions for various rare D meson decays. 
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3 New Physics Analysis 

At this time, there is a wide collection of possible New Physics models leading to FCNC D tran- 
sitions. Among those considered in Ref. 1 are (i) Supersymmetry (SUSY): R-parity conserving, 
R-parity violating, (ii) Extra Degrees of Freedom: Higgs bosons. Gauge bosons, Fermions, Spa- 
tial dimensions, (iii) Strong Dynamics: Extended technicolor, Top-condensation. 

Due to limitations of time (for the talk) and space (for this summary) we restrict most of 
our attention to the case of supersymmetry. However, at the end we also make a few remarks 
on the topic of large extra dimensions. The SUSY discussion divides naturally according to how 
the R-parity Rp is treated, where 

3.1 R-parity conserving SUSY 

R-parity conserving SUSY will contribute to charm FCNC amplitudes via loops. For a penguin- 
like amplitude whose external lines are SM particles, the internal lines can be (i) gluino-squark 
pairs, (ii) charged Higgs-quark pairs or (iii) chargino/neutralino-squark pairs. Case (i) is the 
one considered here as case (ii) will be suppressed for the same CKM reason as in the SM and 
case (iii) is relatively suppressed to (i) because the vertices are weak-interaction rather than 
strong-interaction . 

To calculate R-parity conserving SUSY contributions, we employ the so-called mass inser- 
tion approximation, which is oriented towards phenomenological studies and is also model 
independent. Let us first describe what is actually done and then provide a brief explanation of 
the underlying rationale. 
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Figure 4: A typical contribution to c ^ u FCNC transitions in the MSSM. The cross denotes one mass insertion 
{Si2)xx> and A, A' are helicity labels. 



In this approach, a squark propagator becomes modified by a mass insertion {e.g. the 'x' in 
Fig. ^) that changes the squark flavor. 011 For convenience, one expands the squark propagator 
in powers of the dimensionless quantity i6fj)\y, 
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where i ^ j are generation indices, A, A' denote the chirahty, (Mj")^ are the off-diagonal elements 
of the up- type squark mass matrix and Mq represents the average squark mass. The exchange 
of squarks in loops thus leads to FCNC through diagrams such as the one in Fig. ^. The role of 
experiment is to either detect the predicted (SUSY-induced) FCNC signal or to constrain the 
contributing (^^^Oaa'- 

This topic is actually part of the super-CKM problem. If one works in a basis which diago- 
nalizes the fermion mass matrices, then sfermion mass matrices (and thus sfermion propagators) 
will generally be nondiagonal. As a result, flavor changing processes can occur. One can use 
phenomenology to restrict these FCNC phenomena. The Q = —1/3 sector has yielded fairly 
strong constraints but thus far only D^-D^ mixing has been used to limit the Q = -1-2/3 sector. 
In our analysis, we have taken charm FCNCs to be as large as allowed by the D-mixing upper 
bounds. 



For the decays D Xui'^i discussed earlier in Sect. 2.1, the gluino contributions will occur 



additively relative to those from the SM and so we can write for the Wilson coefficients, 

= C7f^) + Cf . (7) 



To get some feeling for dependence on the ((^i2)aa' parameters, we display the examples 

Cl « (<55'2)ll and (d't^hR , oc (51'2)ll , (8) 

whereas for quark helicities opposite ^ to those in the operators of Eq. (^) , one finds 

cl oc (<5r2)RR and {6^^2)1^^ , C'l cx ((^r2)RR • (9) 

Moreover, the term in which contains ((5j'2)lr experiences the enhancement factor Mg/rric. 

We have numerically studied the effects in c ^ ut^l~ for the range of masses: (I) Mg = Mq = 
250 GeV, (II) Mg = 2Mq = 500 GeV, (III) Mg = Mg = 1000 GeV and (IV) Mg = {1/2) Mq = 
250 GeV. For some modes in D — > the effect of the squark-gluino contributions can be 

large relative to the SM component, both in the total branching ratio and for certain kinematic 
regions of the dilepton mass. The mode p^e^e~ is given in Fig. ^. This figure should 

make clear the importance of measuring the low m^+g- part of the dilepton mass spectrum. 



"^We use the notation C for the associated Wilson coefficients. 
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Figure 5: Dilepton mass distributions for D° — > p'^e^e in the mass insertion approximation of MSSM. The SM 
prediction (solid curve) is provided for reference and the MSSM curves refer to (i) Mg = Mq = 250 GeV, (ii) 
Mg = 2Mq = 500 GeV, (iii) Mg = Adg = 1000 GeV and (iv) Mg = Mq/2 = 250 GeV. 



3.2 R-parity violating SUSY 

The effect of assuming that i?-parity can be violated is to allow additional interactions between 
particles and sparticles. Ignoring bilinear terms which are not relevant to our discussion of 
FCNC effects, we introduce the -R-parity violating (RPV) super-potential of trilinear couplings, 
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where L, Q, E, U and D are the standard chiral super-fields of the MSSM and i,j,k are 
generation indices. The quantities Xijk, X'ijk and are a priori arbitrary couplings which 
total 9 + 27 + 9 = 45 unknown parameters in the theory. 

For our purposes, the presence of RPV means that tree-level amplitudes become possible in 
which a virtual sparticle propagates from one of the trilinear vertices in Eq. (|l^) to another In 
order to avoid significant FCNC signals (which would be in contradiction with current experi- 
mental limits), bounds must be placed on the (unknown) coupling parameters. As experimental 
probes become more sensitive, the bounds become ever tighter. In particular, the FCNC sector 
probed by charm decays involves the {A^^^}. Introducing matrices Ui,, to rotate left-handed 
up-quark fields and right-handed down-quark fields to the mass basis, we obtain for the relevant 
part of the superpotential 

= K,k [-eld\u{ - u^Sl^el - {Sknelfui + ..] , (11) 

where neutrino interactions are not shown and we define 

Some bounds on the {X'ij^} are already available from data on such diverse sources as charged- 
current universality, the ratio FTr-^ei^e /Ftt^^;/^ , the semileptonic decay D —>■ K£u£ , etc.i^ We 
have considered additional experimental implications of the preceding formalism: 

(i) For the decay Tr~^e^e~ , we display the effect of RPV as the dot-dash line in Fig. |l|. 

Here, the effect is proportional to X[^j^ ■ X[2k and we have employed existing limits on these 
couplings. Although the effect on the branching ratio is not large, but the dilepton spectrum 
away from resonance poles is seen to be sensitive to the RPV contributions. This case is not 
optimal because the current experimental limit on i3r£)+_,^+g+g- is well above the dot-dash 
curve. 



(ii) For n^fi'^n , the current experimental limit on Br£)+^^+^+^- actually provides 

the new bound 

A'm • A'm < 0.004 . (13) 



(iii) Another interesting mode is fi^fi . Upon using the bound of Eq. (13) we obtain 



Brf,l + _ < 3.5 X 10-6 ( ^ 



' I 0.02 / 



(14) 



A modest improvement in the existing limit on Brjjo^^j_+^- will yield a new bound on the product 

(iv) Lepton flavor violating processes are allowed by the RPV lagrangian. One example is 
the mode — > e~^n~, for which existing parameter bounds predict 
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An order-of-magnitude improvement in Br 
combination of RPV couplings. 



will provide a new bound on the above 



3.3 Large Extra Dimensions 

For several years, the study of large extra dimensions ('large' means much greater than the 
Planck scale) has been an area of intense study. This approach might hold the solution of the 
hierarchy problem while having verifiable consequences at the TeV scale or less. Regarding the 
subject of rare charm decays, one's reaction might be to ask How could extra dimensions possibly 
affect the decays of ordinary hadrons?. We provide a few examples in the following. 

Suppose the spacetime of our world amounts to a 3 + 1 brane which together with a manifold 
of additional dimensions (the bulk) is part of some higher-dimensional space. A field Q which 
can propagate in a large extra dimension will exhibit a Kaluza-Klein (KK) tower of states {On}) 
detection of which would signal existence of the extra dimension. Given our ignorance regarding 
properties of the bulk or of which fields are allowed to propagate in it, one naturally considers 
a variety of different models. 

Assume, for example, the existence of an extra dimension of scale 1 /R ~ 10^^ eV such that 
the gravitational field (denote it simply as G) alone can propagate in the extra dimension. 
There are then bulk-graviton KK states which couple to matter. In principle there will be 

the FCNC transitions c — > n G„ and since the {Gn} remain undetected, there will be apparent 
missing energy. However this mechanism leads to too small a rate to be observable. 

Another possibility which has been studied is that the scale of the extra dimension is 1/ii ~ 
1 TeV and that SM gauge fields propagate in the bulk. EJ However, precision electroweak data 
constrain the mass of the first gauge KK excitation to be in excess of 4 TeVllj, and hence their 
contributions to rare decays are small li^. 

More elaborate constructions, such as allowing fermion fields to propagate in the five- 
dimensional bulk of the Randall-Sundrum localized-gravity modellll, are currently being actively 
explored. c3 Interesting issues remain and a good deal more study deserves to be done. 



4 Concluding Remarks 

In the Standard Model, FCNC processes are suppressed. Both photon and Z-boson SM vertices 
are flavor diagonal, so tree level diagrams involving virtual propagation of these particles will not 
contribute to FCNC processes. One must instead consider loops. As regards FCNC processes, 



SM loops are largest for kaon and B-meson transitions due mainly to the large t-quark mass and 
also to favorable CKM dependence. The upshot is that studies of FCNC processes for charm 
have lagged behind those of the other flavors. There is of course a basic irony in this - it is 
precisely because SM signals are expected to be so small for charm FCNC that the opportunity 
for evidence of New Physics to emerge becomes enhanced relative to the other flavors. The 
situation begs for experimental input. 
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